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ABSTRACT 

We analyse the evolution of the Intergalactic Medium (IGM) by means of an ex- 
tended set of large box size hydrodynamical simulations which include pre-heating. 
We focus on the properties of the z ~ 2 Lyman-a forest and on the population of 
clusters and groups of galaxies at z = 0. We investigate the distribution of voids in 
the Lyman-a flux and the entropy-temperature relation of galaxy groups, comparing 
the simulation results to recent data from high-resolution quasar spectra and from 
X-ray observations. Pre-heating is included through a simple phenomenological pre- 
scription, in which at z — 4 the entropy of all gas particles, whose overdensity exceeds 
a threshold value <5h is increased to a minimum value Kq. While the entropy level 
observed in the central regions of galaxy groups requires a fairly strong pre-heating, 
with Ka > 100 keV cm 2 , the void statistics of the Lyman-a forest impose that this 
pre-heating should take place only in relatively high-density regions, <5h~ 30, in order 
not to destroy the cold filaments that give rise to the forest. We conclude that any 
injection of non-gravitational energy in the diffuse baryons should avoid low-density 
regions at high redshift and/or take place at relatively low redshift zlS 1. 
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1 INTRODUCTION 

Observations in the X-ray band of the hot intra-cluster and 
intra-group medium, along with observations of the absorp- 
tion features in the spectra of distant quasars (QSOs) from 
the intervening intergalactic medium (IGM) offer powerful 
means of tracing the evolution of diffuse cosmic baryons in 
different regimes. While X-ray observations of galaxy sys- 
tems trace the high-density bar yons in the low-redshift z^, 1 
Universe (e.g. IVoit e t al. 2005), data on the Lyman-a for- 
est convey informatio n on baryons a round the mean cosmic 
density at z^ 2 (e.g. iMei ksin 2007). As such, these differ- 
ent observational techniques complement each other in the 
reconstruction of the cosmic cycle of baryons. 

One of the standard results from the X-ray observa- 
tions of groups and clusters of galaxies is that the hot gas 
in the central regions of low-temperature systems has an 
entropy level higher than predict ed by the gravitationa l 
process of accretion shocks (e.g. iTozzi fc Normanl 1200 ll ). 
The commonly accepted explanation for this is that radia- 
tive cooling and heating from some feedback energy source 
should be the main mechanisms responsible for setting the 
hot intra-cluster medium on a relatively high adiabat, while 
preventing overcooling from co nverting into star s a unrealis- 
tically large fraction of baryons (|Voit et al.ll2005l ) . In order to 



test these predictions, hydrodynamical simulations includ- 
ing some form of e xtra gas heating hav e been analysed by a 
number of authors (|Borgani et al.l 2008, for a recent review). 
The general result from these simulations is that injecting 
~ 1 keV per gas particle is effective in increasing the central 
gas entropy to a level consistent with the observed one, while 
the question remains open as to what astrophysical mech- 
anism should be responsible for this high-redshift heating. 
Clearly, any diffuse IGM heating at high redshift is expected 
to leave an imprint on its observational properties. For in- 
stance, increasing the gas temperature inside the filaments 
permeated by the IGM should significantly alt er the statis- 
tical p roperties of the Lyman-a forest. Indeed, [Shang et al.l 
(2007) recently suggested that the pre-heating required to 
reproduce X-ray cluster observations should leave the im- 
print of bubbles of ionised gas at high redshift. These authors 
analysed a sample of SDSS QSO spectra searching for large 
voids in the Lyman-a forest associated to these bubbles. 
From the lack of detection of these voids, they concluded 
that any pre-heating should only involve IGM at a density 
higher than the mean one as traced by the Lyman-a forest. 

In this Letter we present a combined analysis of the 
void statistics in the Lyman-a forest at z = 2.2 and of the 
entropy level of the hot gas within galaxy groups at 2 = 0, 
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Figure 1. The different panels show the gas overdensity in a 
slice of 8 comoving /i _1 Mpc for the reference simulation, with no 
pre— heating (left panels) and for the (300,10) simulation (entropy 
floor of 300 keV cm 2 for (Sgas > 10, right panels). Top and bottom 
panels refer to the snapshots at z = 2.2 and at z = 0, respectively. 



from a set of pre-heated cosmological hydrodynamical sim- 
ulations within a fairly large cosmological box. The analysis 
of the void statistics in the Lyman-a forest will be carried 
out b y using the same procedure presented by IViel et al.l 
(2008), where a description of the reference observational 
data set is also provided. As for the entropy of the intra- 
group gas, we will compare simulation predi ctions to the re- 
cently published results bv lSun et al.l f2008), who analysed 
CHANDRA data for an extended set of nearby groups. By 
changing both the strength of the pre-heating and the min- 
imum overdensity at which such a pre-heating takes place, 
we aim at better exploiting the complementary information 
placed by high-density hot gas in nearby groups/clusters and 
by the lower-density IGM associated to the high-redshift 
Lyman- a forest. 



2 HYDRODYNAMICAL SIMULATIONS 

We use simulations car ried out with the Tree-SPH GADGET- 
2 code (ISpringell fcoOS). These simulations include radiative 
coolin g and heating from a uniform evolving UV background 
(see iBolton et all ( 20051 ) for details), for a primordial mix 
of hydrogen and helium. Since no self-consistent model for 
chemical enrichment is treated in these simulations, we de- 
cided not to include the effect of metals in the computation 
of the cooling function. The star formation criterion sim- 
ply converts in collisionless stars all the gas particles whose 
temperature falls below 10 K and whose densit y contrast 
is larg er than 1000. More details can be found in 1 Viel et al.l 
( 2004). The cosmological model assumed for our simulations 
corresponds to a ACDM Universe with f2 m = 0.26, Qa = 



and erg = 0.85, in agreement with the latest results of cosmo- 
logical parameter s derived from CMB and l arge-scale struc - 
ture observables (jLeseourgues et al"1l2007l : lKomatsul [2008). 
We have used 2 x 400 3 dark matter and gas particles in a 
120 ft" 1 comoving Mpc box for our purposes. The gravita- 
tional softening was set to 15 ft" 1 kpc in comoving units 
for all the particl es. Even if these s imulations barely resolve 
the Jeans length, IViel et all (|2008l ) showed that the statis- 
tics of voids have numerically converged once the simulated 
spectra have been smoothed over a scale of 1 ft Mpc co- 
moving. With this mass resol ution, the smalle st groups in 
the observational data set by ISun et al.1 Hool) which have 
are resolved with more than 5000 par- 
ticles. 

Non-gravitational heating is included following a stan- 
dard purely phenomenological model of pre-heating, based 
on imposing a minimum entropy floor, Kn to all the gas 
particles that, at a given heating redshift Zh, have density 
contrast above a given threshold S^. We use the definition of 
entropy which is standard in the study of the intra-cluster 
medium, namely, K — n e /T, where n e is the electron num- 
ber density and T the gas temperature. Similar schemes 
to study the effect of pre-heating to the X-ray properties 
of galaxy clusters have been used by several authors (e.g. 
i Bialek et"aL 2001 ; Muanwong et al. 20021; iTornatore etall 



20031 : iBoreani et al]|2005l :TY< 



0.74, Q b = 0.0463, 



0.95, and H = 72 km s" 1 Mpc" 



fc Brvanll2007h . The pn> 
cedure adopted to pre-heat the gas is the following: at z^ = 4 
we select from the reference run, which does not include any 
pre-heating, all the gas particles having 5 > 5h and en- 
tropy K < Kfi- The internal energy of these particles is 
then increased so that their entropy match the floor value. 
In this selection, we always exclude cold and dense star par- 
ticles, with 8 > 500 and T < 3 x 10 4 K, whose low ini- 
tial entropy would require an extremely strong and ad-hoc 
heating. Such cold and dense gas is contained wit hin galaxy- 
sized halos, where star formation takes place (e.g., lKav et al.l 
2000). Therefore, imposing the above limits on density and 
temperature of the gas to be pre-heated amounts to protect 
the cold gas content of galaxies, which should be rather af- 
fected by some local feedback mechanims (e.g., associated 
to SN explosions). We refer to iVoitl (|2005l ) for an overview 
of the role of pre-heating in modelling the thermodynamics 
of the intra- group medium. 

We choose four different values of the entropy floors: 
Kd — 10, 50, 100, 300keVcm 2 ; and 3 overdensity thresholds: 
<5h = —1,3,10. Only for the case with Ka — 300keVcm 2 
we also consider the case with a higher density threshold for 
heating, <5h = 30. While the lowest <5h corresponds to heating 
all the gas, the highest value is comparable to that expected 
at the boundary of a virialized structure. In this way, we run 
a total of 14 simulations (13 different pre-heated runs and 
the default run with no pre-heating). In the following, we 
will show results only for a subset of the whole simulation 
set. The different simulations will be labelled by the tuple 
of values of entropy floor (in keVcm 2 ) and of overdensity 
threshold for heating. The reference run will be indicated 
with the label REF. All the simulations are analysed at z — 
2.2 to compute the Lyman-a forest statistics and at z — 
to compute the entropy-temperature relation of the galaxy 
groups. 

The sample of high -resolution, high s ignal-to-noise QSO 
spectra is described in iKim et al.1 l|2007h and has been al- 
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Figure 2. The ratio between Lyman-a flux voids larger than Rfi ux (in comoving h _1 Mpc) of pre— heated and the no— heating (REF) 
case. In the left panel we show the effect of a different gas density threshold having fixed the entropy floor for the models (300,-1), 
(300,3), (300,10) and (300,30). In the right panel we show the effect of a different value for the entropy floor having fixed the density 
threshold for the models (10,10), (50,10), (100,10), (300,10). The shaded area in both panels could be considered as the 1 — a estimate 
of the observational statistical error (see text). 



ready used to compute seve ral flux statistics (jBolton et al.l 
12008 : IViel et alJl2004L l2008h . Many systematic effects have 
been carefully addressed including metal lines contamina- 
tion and continuum fitting errors. The median redshift of 
the sample is at z = 2.2. From this data set we extracted 
the cumulative distribution of voids in the transmitted flux, 
i.e. connected regions with flux above the mean level, which 
at 2 ~ 2 appear to trace reasonably well underdense regions. 

As for the CHANDRA sa mple of ga l axy g roups, it con- 
tains 40 objects selected by ISun et all l|2008f) so that gas 
properties can be derived out to r25oo for all of them (Ya 
is the radius encompassing an aver age density of A x p c , 
being p c the cosmic critical density). ISun et"aL (|200&f ) anal- 
ysed this sample to derive the scaling relation between en- 
tropy and temperature at different overdensities and found 
an excess of entrop y at rspp, w i th re spect to the baseline 
value calibrated by IVoit et ail (|2005h from non-radiative 
hydrodynamical simulations of clusters, although this ex- 
cess is smaller than suggested by previous analyses (e.g. 
Ponman et all l2003l : iPratt fc Arnaudll2005l : iPiffaretti etatl 



20051 ). Moreover, they also found that the excess entropy is 



larger at r2soo, thus confirming that any non-gravitational 
process has a larger effect in the central regions of groups. 



3 RESULTS 

In Figure [T] we show a projected slice whose thickness is 8 
comoving /i~ 1 Mpc of the gas distribution at z = 2.2 and 
2 = 0, for the default run and for the (300,10) run. Even by 
adding this strong pre-heating, at z = 2.2 the skeleton of the 
cosmic web at densities around the mean is still preserved. 
The main differences arise in dense structures, whose filling 
factor is small. At z — the effect of pre-heating is quite 
visible as well: the mildly non-linear cosmic web evolves and 
gives rise to clusters of galaxies and galaxy groups. These 
very non-linear structures tend to be punier and smoother 
compared to the reference case, with a suppression of the 
number of small halos that trace the filaments. 



3.1 Void statistics in the Lyman-a forest 

From the snapshots at z — 2.2 we extract a mock set of 
1000 Lyman-a QSO spectra in random directions. Then, we 
smooth the spectra over a scale of 1 comoving /i -1 Mpc, to 
be less sensitive to small structures at and around the Jeans 
length that might not be properly resolved by our simu- 
lations. We then compute the number of voids in the flux 
distribution, having size larger than a given value and com- 
pa re the pre-heated runs with the default one. As shown 
by IViel et ail (|200ST ). the reference run provides results in 
excellent agreement with the data. Results are presented in 
Figure [2] where the shaded area indicates the uncertainty 
in the observed mean flux level, which sets the criterion for 
the selection of voids. Among all the possible uncertainties 
in the astrophysical and cosmological parameters , this error 
has t he largest effect on the void statistics (see IViel et ail 
120081 for more details) . Therefore, we take it as a rough es- 
timate of the error. We remind that the total error budget 
must also take into account the contributions from all the 
other parameters and it is larger than the one shown here. 
For example, a ~ 3 times higher temperature of the IGM at 
the mean density would boost the number of 30 comoving 
/i -1 Mpc voids up to a factor 1.5, bringing some of the mod- 
els in better agreement with the observations. Basically, the 
error bars represented by the shaded area do not take into 
account different thermal histories for the low density IGM 
and/or different cosmological scenarios (warm dark matter 
or extra power at intermediate scal es) that are more exten- 
sively discussed in IViel et all (|2008l ). 

In the left panel, we address the role of a different over- 
density threshold for heating, keeping fixed the entropy floor 
at 300 keV cm 2 . Clearly, heating up the whole IGM (i.e. 
<Sh = — 1) substantially increases the number of voids since 
the neutral fraction will be reduced by the very large amount 
of heating. If we increase the overdensity threshold to Sh — 3 
and 10, then the void fraction is in better agreement with the 
default case, as expected, but there is a tendency to under- 
predict the number of large void regions. This result might 
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seem counter-intuitive, since a heating of the IGM should 
reduce the hydrogen neutral fraction. However, we found 
that the gas in the relatively low density environments of 
the pre-heated runs is denser compared to the default run 
and thereby carries with it a larger neutral hydrogen frac- 
tion. The opposite trend takes place in overdense regions 
where the default simulation is denser than the correspond- 
ing pre-heated ones. At low redshifts the volume filling fac- 
tor of underdense (overdense) regions gets larger (smaller) 
in all the models but the overall effect is a reduction of the 
size of large voids compared to the REF case since in the 
pre-heated runs the voids are less empty. The (300,30) is 
closer to the default case than the (300,10) owing to the de- 
crease of the volume filling factor of the heated regions with 
increasing S c . 

In the right panel of Figure [5] we show the effect of 
changing the entropy floor, while keeping the overdensity 
threshold for heating fixed at 5h = 10. Increasing the entropy 
floor has the effect of lowering the fraction of large voids: at 
z — 4 the heating of dense gas particles at S > 10 makes 
them leave quickly their halos and reach the low density 
IGM: these particles have usually a larger fraction of neutral 
hydrogen that can cause absorptions in the mock spectra 
that are extracted at z = 2.2. 



3.2 The entropy— temperature relation of groups 

At z = we identify galaxy groups in our simulations 
by running a friends-of-friends algorithm with b — 0.15 
for the linking length in units of the mean dark matter 
(DM) interparticle separation. The centre of each group is 
then identified with the position of the DM particle having 
the minimum va lue of the gravitational potential. Following 
ISun et all |2008t ). we compute for each group the tempera- 
ture Ta within ta (for A = 500 and 2500) by excluding the 
core region within 0.15 rsoo. Then we measure the profiles of 
the electron number density, n e (r), and of the gas temper- 
ature, T(r). The values of entropy Ka at ta are then com- 
puted as Ka = Ta /i\ a ■ -X-r ay tempera t ures a re computed 
following the prescription by IVikhlininl (|2006h , which rep- 
resents an extension to low-temperature systems (Tx < 3 
keV, relevant for our analysis) of t he spectroscopic-like te m- 
perature, originally introduced bv lMazzotta et ail (|2004l ). 

In Figure [3] we compare the relation between entropy 
and temperatu re for our simula ted groups to the observa- 
tional data by ISun et alj (|2008l ). for four of our simulated 
boxes. As for the case with no pre-heating (upper left panel), 
a satisfactory agreement with data is obtained at rsoo, while 
simulated groups have a too low level of entropy at r2soo- 
This resu lt is in line with the comparison performed by 
ISun et all l|2008h between their observational results and 
th e best-fit rel a tion f rom the hydrodynamical simulations 
bv lNagai etal] (|2007T l. These simulations include radiative 
cooling, star formation and a rather inefficient form of feed- 
back, and therefore are expected to provide similar results 
to the reference (non pre-heated) run. Clearly, the enntropy 
level at r2soo is not the only problem suffered by a radiative 
run without extra heating. Indeed, in this simulation over- 
cooling causes about 40-50 per cent of the baryons within 
rsoo to be converted in stars, with a decreasing trend with 
the system temperature. This fraction, which is a lower limit 



owin g to the resolution dependence of the coo ling efficiency 
(e.g. iBalogh et ail l200ll : iBorgani et ail |2006h. is m excess 
with respect to observational estimates fe.g. iGonzalez et al.l 
[2007|) . Therefore, whi le cooling plays the r ole of establishing 
the level of entropy (|Voit fc Brvanl l200ll ). a form of non- 
gravitational heating is required to regulate the amount of 
lower-entropy gas which is destined to cool and form stars. 

As for heating with Ka — 100 keV cm 2 (upper right 
panel of Fig. [3J), it has a rather small effect on the gas en- 
tropy, while it does have a significant effect on the fraction 
of stellar mass, which drops to 10-20 per cent. This result 
implies that, with this level of Kg., radiative cooling is still 
the main responsible for setting the entropy level, while ex- 
tra heating regulate the amount of cooled gas. A further in- 
crease of iffl is then required to alleviate the tension between 
the observed and the measured levels of entropy in galaxy 
groups. Indeed, using Ka = 300 keV cm 2 (bottom panels 
of Fig. |3J) brings the entropy level in the simulated groups 
to better match the observed one. However, as discussed 
in the previous section imposing an entropy floor generates 
too large voids in the Lyman-a forest, an effect that can be 
compensated by increasing the heating overdensity thresh- 
old, Sh- However, increasing the latter from Sh = 10 (bottom 
left panel) to 30 (bottom right panel) induces a slight but 
sizable decrease of Ka, especially for A = 2500, as a con- 
sequence of the smaller number of gas particles heated at 
z = 4. Therefore, while heating at relatively high overden- 
sity is required by the void statistics of the Lyman-a forest, 
it goes in the wrong direction to reproduce the thermody- 
namical properties of intra-group medium at small radii. 
This demonstrate how effective is combining observations of 
the high-2 IGM and low-z intra-group medium to constrain 
the thermal history of the cosmic baryons. 



4 CONCLUSIONS 

We presented an analysis of cosmological hydrodynamical 
simulations, including radiative cooling and pre-heating, 
aimed at characterising both the properties of the high- 
redshift Lyman-a forest and the thermodynamical prop- 
erties of the diffuse gas within nearby galaxy groups. We 
use a simple phenomenological recipe for pre-heating, in 
which at z — 4 all the gas particles lying at an overdensity 
above 8^ are brought to a minimum entropy level of Ka- At 
z = 2.2 mock QSO spectra were extracted and their prop- 
erties compare d to the observe d distribution of voids in the 
Lyman-a flux (|Viel et al.ll2008l ). while at z — the entropy- 
temperature relation of galaxy groups is compare d with re- 
cent results from X-ray CHANDRA observations l|Sun et al.l 
2008). The main results of our analysis can be summarised as 
follows: i) pre-heating all the gas, irrespective of its density, 
produces voids in the Lyman-a forest, which are too large if 
compared to observations; ii) imposing a modest overden- 
sity threshold for heating, 5^ — 3, suddenly reduces the size 
of the voids to values which are even too small; Hi) fur- 
ther increasing 8^ makes the void sizes approaching those of 
the non pre-heated run, owing to the progressively smaller 
amount of heated gas which ends up in the Lyman-a for- 
est; iv) the entropy level within galaxy groups at rsoo can 
be already matched in non pre-heated simulations, which 
however fail at accounting for this level at r2soo; v) a fairly 
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Figure 3. The relation between entropy and temperatur e for simulated groups at z = (small circles) and for the CHANDRA 
observational of nearby groups analysed by llSun et al.ll200SI) (points with error bars). The four panels correspond to the REF run (with 
no pre— heating; top left), to the (100,10) run (top right), to the (300,10) run (bottom left) and to the (300,30) run (bottom right). The 
top and the bottom parts of each panel refer to the relation for A = 500 and A = 2500, respectively. 



strong entropy injection, with Kg, > 100 keV cm 2 , is re- 
quired to match the entropy-temperature relation of poor 
systems at r25oo- 



of heating from astrophysical sources, such as supernovae, 
active galactic nuclei and DM annihilation. 



The need of reproducing the entropy structure of the 
low-redshift intra-group medium and the void statistics of 
the high-redshift Lyman-a forest leads us to conclude that 
for a mechanism of non-gravitational heating to work, it 
must provide a fairly large amount of extra entropy at rela- 
tively high overdensities, comparable to those characteristic 
of virialized halos. Our conclus ion, in agreement with the 
argument bv lShang et"ail (|2007h . is that the amount of pre- 
heating required by low-redshift galaxy clusters and groups 
should avoid low density regions for it not to produce too 
large voids in the Lyman-a forest or should act at lower 
redshift, 1. Admittedly, our scheme of non-gravitational 
heating is an oversimplified one. For this reason, our aim is 
not to look for the best values of Kg and <5h, which are able 
to fit at the same time the observed intra-group entropy 
and void statistics of the Lyman-a forest. Rather, the main 
goal of our analysis is to provide an indication of the general 
properties that a plausible mechanism of non-gravitational 
heating should have. We defer to a forthcoming analysis a 
detailed study of the combined constraints of the low and 
high-z IGM properties using physically motivated models 
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